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1^ ! Abstract 
> , 

I We show that Q-balls naturally exist in the Minimal Supersymmetric Standard 

Model (MSSM) with soft SUSY breaking terms of the minimal N=l SUGRA 
type. These are associated with the F- and D-flat directions of the scalar po- 
, tential once radiative corrections are taken into account. We consider two dis- 

^ I tinct cases, corresponding to the " HuL" (slepton) direction with L-balls and the 

"n'^d'^d'^" and '''' u^u'^iFe"^ (squark) directions with B-balls. The L-ball always 
has a small charge, typically of the order of 1000, whilst the B-ball can have 



in 



O , an arbitrarily large charge, which, when created cosmologically by the collapse 

of an unstable AfHeck-Dine condensate, is likely to be greater than 10^^. The 



X 
J3 



B-balls typically decay at temperatures less than that of the electroweak phase 
transition, leading to a novel version of Affleck-Dine baryogenesis, in which the B 
asymmetry comes from Q-ball decay rather than condensate decay. This mecha- 
nism can work even in the presence of additional L violating interactions or B—L 
conservation, which would rule out conventional Affleck-Dine baryogenesis. 
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1 Introduction 



Flat directions in the scalar potential of the Minimal Supersymmetric Standard Model 
(MSSM) [m have long been understood to have potentially important consequences 
for cosmology ^. In particular, a very natural mechanism for the generation of 
the observed baryon asymmetry in SUSY models is the decay of a charged scalar field 
condensate (the Affleck-Dine (AD) mechanism [0). Developments in the cosmology of 
N=l supergravity models, in particular the generation of order H contributions to the 
soft SUSY breaking parameters of scalar fields in the early Universe , have sparked 
a revival of interest in flat directions and the AD mechanism 0. This mechanism is 
particularly important as the other natural possibility for baryogenesis in the MSSM, 
namely electroweak baryogenesis, appears to succeed for only a very limited range of 
MSSM parameters 0. A second idea has recently been reconsidered in the context of 
SUSY models with flat directions in their scalar potentials, namely Q-balls [7-11]. A 
Q-ball is a stable, charge Q soliton in a scalar fleld theory with a spontaneously broken 
global U{1) symmetry [^]. The connection between Q-balls and SUSY models lies in 
the natural existence of scalars carrying a global U{1) charge, corresponding to B or 
L for the case of squarks and sleptons. Thus, in principle, one could have B- or L-balls 
in the MSSM. Whether this happens or not depends on whether the scalar potential 
of the squarks and sleptons can support a Q-ball solution. If it can, then the next 
questions are whether Q-balls can be created in the early Universe and whether the 
resulting Q-balls could be sufficiently long-lived to have interesting consequences for 
cosmology. In reference |]^ it was shown, in the context of a SUSY model with gauge 
mediated SUSY breaking, that the scalar potential along a flat direction of the MSSM 
can support a Q-ball solution. In addition, it was suggested that Q-balls with a very 
large charge (and so possibly stable) could be formed by the growth of perturbations 
in a charged condensate of the kind associated with the AD mechanism. 

In this letter we will reconsider these ideas in the context of the "classical" MSSM, 
i.e. the MSSM with soft SUSY breaking terms due to SUSY breaking in a gravita- 
tionally coupled hidden sector ||l|. We will flrst show that the scalar potential along 
the flat " H^L" (slepton) and "u^d'^d'^" and "u'^u'^d'^e'^" (squark) directions can indeed 
support a Q-ball solution once radiative corrections are taken into account. The form 
of the potential and the associated Q-balls will be seen to be quite different in the 
slepton and squark cases. We will then discuss the formation of Q-balls in the early 
Universe in the context of the simplest cosmological scenario, in which the AD fleld 
begins coherently oscillating during the inflaton matter dominated era following an 
initial period of inflation. Finally, we will discuss the lifetime of the Q-balls and their 
possible consequences for cosmology, in particular their role in baryogenesis. 



2 Q-balls and radiatively corrected flat directions 
in the MSSM 



2.1 Flat directions in the MSSM 

For a globally U{1) symmetric complex scalar field 0, the condition for the scalar 
potential f/(0) to be able to support a Q-ball solution is that ?7(0)/|0p should have a 
global minimum at a non-zero value of\(f)\ [0 . In this it is assumed that (p is defined such 
that = corresponds to the vacuum in which the Q-ball is formed, which in practice 
corresponds to our vacuum with non-zero Higgs expectation values and zero squark and 
slepton expectation values. We will be interested in whether Q-balls can be formed 
along flat directions of the MSSM. The flat directions can be defined by the scalar 
field operators which have non-zero expectation values along the directions in question. 
These have been classified generally for the case of the MSSM in reference . The most 
interesting fiat directions from the point of view of conventional AD baryogenesis are 
those with non-zero B-L and unbroken R-parity, which can give a baryon asymmetry 
once anomalous B+L violation is taken into account without introducing dangerous 



renormalizable B and L violating operators in the MSSM [||, |12[. These fiat directions 
correspond to the " HuL" direction, along which the magnitude of the expectation value 
of and ul are equal, and the u'^d'^d'^" direction, along which Uc, dc and d'^ are non- 
zero, where the squarks have different colour indices and where dc and d'^ correspond 
to orthogonal combinations of down squark generations. (There are also the "d'^QL" 
and "e^LL" directions. However, since we expect these to be phenomenologically 
similar to the "u'^d'^d'^" direction, we will not consider them seperately here). R-parity 
conservation allows the d=4 non-renormalizible superpotential term (if^L)^ and the 
d=6 term {u'^d'^d^Y. In addition, we will consider the d=4 B-L conserving u^u'^d'^e'^ 
direction (there is also the phenomenologically similar QQQL direction). Although 
this could not produce a B asymmetry via AD condensate decay once anomalous B + L 
violation is taken into account, we will see that it can generate a baryon asymmetry 
via Q-ball decays occuring after the electroweak phase transition. 

For sufficiently large |0| along the D-fiat direction the scalar potential is of the 
form 



,2 , A^l^pC'^-i) , ,AaA0^ 



IVlp 

where m| is the soft SUSY breaking mass squared term and Ax the A-term, and we 
take the non-renormalizible terms to be suppressed by the Planck scale Mp. Thus, if 
has no |0| dependence, then f/(0)/|0p will not have a global minimum at non-zero |0| 
and so no stable Q-ball will exist. In order to have such a |0| dependence, Kusenko and 
Shaposhnikov considered the possibility that m| is generated by gauge mediated 



SUSY breaking. In this case, for \(f)\ larger than the mass of the messenger quarks, mg 
becomes proportional to and the potential becomes essentially completely flat up 
to where the non-renormalizible terms become important. However, gauge mediated 
SUSY breaking is not the most commonly considered form of SUSY breaking. The 
most common SUSY breaking mechanism is that where soft SUSY breaking terms are 
generated by SUSY breaking in a hidden sector, not coupled directly to the observable 
MSSM fields in the Kahler potential of N=l supergravity This results in soft SUSY 
breaking terms which are constant at tree level. However, once radiative corrections 
are taken into account, they become scale dependent. Indeed, it is widely believed that 
radiative corrections to the soft SUSY breaking mass squared term of the Higgs scalar 
giving masses to the up-type quarks, Hu-, when summed from the GUT or Planck scale 
using the renormalization group (RG), are responsible for driving its mass squared 
negative at renormalization scales of the order of the weak scale, so breaking the 
electroweak symmetry |13[. It is well-known that the effect of radiative corrections 
on the effective potential as a function of |0| can be summed by replacing the tree- level 
masses and couplings by their value run by the RG equations to a renormalization scale 
/i ~ 101 Thus, for sufficiently small values of |0|, the scalar potential along the 



D-flat direction will have the form U{(j)) ~ m|(|0|)|(;/)p. A necessary condition for this 
potential to be able to support a Q-ball solution is that m|(|0|) decreases over at least 
some range of \(f)\ as |0| increases from zero. 



2.2 HuL direction 

In general we will be considering models with several scalar fields. To correctly describe 
the Q-balls, which correspond to L-balls in this case, we would have to solve the coupled 
equations of motion of the scalars to find a Q-ball solution. However, in practice, we 
will be able to describe the nature of the Q-balls reasonably accurately by considering 
a single scalar degree of freedom. This can be understood as follows. Along the D-fiat 
direction, for large values of the scalar field compared with the mass scale of the soft 
SUSY breaking terms (which are typically of the order of 100 GeV), the degree of 
freedom in the direction orthogonal to the D-fiat direction will be much more massive 
than 100 GeV and so will play no role in the dynamics of the Q-ball. Thus we need 
only consider the scalar field with non-zero value along the D-fiat direction. However, 
at sufficiently small values of the scalar field, the orthogonal scalar mass becomes of 
the same order as the soft SUSY breaking masses and so all the scalar fields can play 
a role in the dynamics of the Q-ball. However, in this limit, all the fields will typically 
have similar masses and couplings and so we can simply consider one of the scalars 
to be decoupled from all the others and calculate the Q-ball dynamics for this scalar. 
This will give the correct order of magnitude for the properties of the Q-ball in this 
regime. 



Let us first consider the potential at large enough values of \(f)\ that the HuL D-fiat 
direction describes the minimum of the potential. In this case the scalar field along 
the minimum of the potential corresponds to a linear combination of the H° and z/^ 
fields, (j) = -^{H^ + ^l), with mass m| = |(yu|^ + + where fin corresponds 
to the SUSY /i-term in the superpotential, iinHuHd- In general, the renormalization 
group equations for the soft SUSY breaking mass terms will have the form 

din ' 

= a(m\ - ^^Ml , (2) 

where the represent the soft SUSY breaking scalar mass terms as well as the A- 
terms, and are the gaugino masses. The full RG equations are given in reference 
[0]. The only important terms are those associated with the top quark Yukawa 
coupling. These appear in the RG equation for mu^ but not for mi^. In the absence 
of a large the masses increase monotonically with decreasing /i. However, with a 
large ctj, for sufficiently small /z, the masses can start to decrease with decreasing [i. 

The form of the solutions depends on the initial values of the parameters of the 
SUSY breaking terms at the initial large scale, which we will take in the following to 
be the GUT scale Mx- However, two features are important for the existence and 
cosmology of L-balls. First, in general, the value of m| becomes negative for scales 
typically smaller than 10*^ GeV or so. Second, depending on the parameters, there can 
be a "hill" in the plot of m| versus |0|, such that m| starts decreasing with increasing 
1 01 for sufficiently large values of |0|. 

The running of the mass in the i/^jL-direction as obtained from solving the coupled 
one-loop RG equations is illustrated in Fig. la for different initial values of the SUSY 
parameters, chosen so that they give radiative breaking of electroweak symmetry at 
the scale /i ~ M^y (we took as{Mw) = 0.11 and the top Yukawa ht{Mw) = 1.05; no 
threshold corrections were augmented). We have not made a systematic exploration of 
the parameter space, but it appears that typically the "hill" in the plot of m| versus 
101 arises for mo smaller than the gaugino masses or for A{Mx) < 0. Examples of 
both can be seen in Fig. la. 

The effect of the negative value of m| at small enough |0| is to generate a minimum 
of [/(|0|)/|0p, typically at |0o| ~ 1 TeV. To see this, we first note that the definition 
of the L-ball field will not correspond to the D-fiat direction once H° ^ Tni/ g. The 
vl expectation value vanishes at the minimum of the potential as a function of if° 
in this case. In this regime the L-ball field may be taken to be vl. This will have a 
positive mass squared, corresponding to our L-conserving vacuum. Thus we see that 
the effective L-ball field will be characterized by a positive mass squared for |0| ^ 1 TeV 
and a negative mass squared for |0| ^ 1 TeV. The rate by which the positive mass 
switches on is given by the rate by which the Higgs VEV develops. This mass turns 
positive over a very narrow range in 0. This is illustrated in Fig. lb. We need not 
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Figure 1: a) The running of the mass in the if^L direction (lower set of curves) for 
the susy parameters (a) mo = 0, A{Mx) — 0, P^h — 2.55, (b) mo = 0, A{Mx) — 
—3, hh — 2.4 and (c) mo = 1, A{Mx) — —1, /jtH — 2.5, all in units of the common 
gaugino mass; for comparision, the running of the physical Higgs mass is also shown 
(upper set of curves), b) The Q-ball potential for the case (a). 

be concerned with the details of the potential; the L-ball properties will be essentially 
determined by this rapid change of the sign of the effective mass squared of the L-ball 
field, which is responsible for f/(|0|)/|(;/)p having a minimum at |0o| ^ 1 TcV. 

The potential importance of the "hill" feature in Fig. lb is that if a condensate in 
the early Universe starts to oscillate coherently with its initial amplitude in the range 
of 101 where m| is decreasing with increasing |0|, then, as we will show later, there 
could, in principle, be a negative pressure associated with the condensate, which would 
lead to the exponential growth of perturbations in the condensate field and ultimately 
to the formation of L-balls. 

2.3 u'^d^d^ direction 

In this case the D-flat direction essentially describes the B-ball field for all values of 
|0|. The corresponding mass term will have the form m| ~ ajm^ ^ + Am^ j + 7im^/ ., 
where the sum is over different generations i and where d and d' correspond to different 
colours and orthogonal combinations of down squark generations, in order to have a 
non-zero value for the operator e^^'^u^djd'^' , where i, j and k are colour indices. The 
values of Pi and will be determined by the initial conditions following inflation, 
which will determine the direction in colour and flavour space of the fiat direction. On 
average, in the sense of allowing equal contributions from all squark species, we expect 
that the effective mass will be of the form m| ~ i(m^c+2m^c+6m^c) (where we assume 



a common mass for first and second generation up squarks and for all down squarks) , 
although it is quite possible to have, for example, a much smaller contribution from 
mfc (which can decrease with decreasing |0| for small enough |0|) for particular initial 
conditions. We generally expect m| to decrease with increasing |0| without any hill 
feature or change of sign, in complete contrast with the case of the H^L direction. 
This leads to a correction to the mass which we can roughly model as a logarithmic 
correction 

rnl^ml(l + K\og(^^^ , (3) 

where < and Mx is a large mass scale at which m| is defined to have the value 
m^. Perturbatively, the value of K due to SUSY breaking gaugino masses will be of 
the form 

K^- V (4) 

a, gauginos S 

which, taking ^0.1 for SU{'i)c and summing over the gauginos gaining a 
dependent mass from mixing with the u'^ and squarks, will typically give in 
the range 0.01 to 0.1. The minimum of the potential will then be determined by the 
non-renormalizible terms in the potential. The form of the non-renormalizible terms 
will depend on the flat direction in question. Since we are considering R-parity to be 
unbroken, the lowest-order non-renormalizible superpotential term lifting the scalar 
potential along the u^d'^d^ direction will be of the form -^{ifd^d'^Y . Thus the full 
scalar potential in this direction will have the form 

U{4>) ^ ml [l + Klog + ■ (5) 

In what follows we will always take A ~ 1. Thus the value of |0| at the minimum of 
[/(|0|)/|0|2is givenby|0|of5i (|K|m|M|;)i/8 5 x lO^^I/Tl^/^ GeV for ~ 100 GeV. 

2.4 u'^u'^d'^e'^ direction 

This will be similar to the u'^d'^d'^ direction, except that the non-renormalizable terms 
will correspond to the d=4 superpotential term -^u'^u'^d'^e'^. Thus the scalar potential 
in this case will have the form 

U{ct>) « ml (l + i^log (^Ig)) 101^ + ^l^r ■ (6) 

The value of |0| at the minimum of C/(|0|)/|0p is given by |0|o ~ (|-fr|m|M|,J^/^ 
3 X 10i°|X|V4 GeV for ms f« 100 GeV. 



3 Q-ball solutions 



3.1 Q-ball equation of motion 



From the point of view of cosmology and phenomenology, the important quantites are 
the energy and radius of the Q-ball as a function of its charge Q. The Q-ball solution 
is of the form 

(7) 



V2 



From now on will refer to 0(r), which is real and positive. The energy and charge 
of the Q-ball are then given by 



E 



2 V dr 
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(8) 



and 



Q 



The equation of motion for a Q-ball of a fixed value of uj is given by 



2 / 
+ -0 
r 



(9f/(0) 



(9) 



(10) 



where 0' = d(j)/dr. We require a solution such that 0(0) = 0'(O) = and ^ as 
r oo. This corresponds to a tunnelling solution for the potential —U{(f)) M, where 



U{<j>) = f/(0) - y02 . 



'Ill 



In practice we vary 0(0) with these boundary conditions until the correct form of 
solution is obtained for a given u. The energy and charge of the solution are then 
calculated using the above expressions. 



3.2 H^L direction 

The scalar potential which describes the L-ball solution for the HuL direction will be 
that due to the rapid change in the sign of the mass squared term of the effective 
L-ball field at |0| < 1 TeV. We will model this by 

2 

771 

f/(0)^_(2e-^_l)02^ (12) 

where s ~ ITeV^^. It is straightforward to show (and confirmed by numerical solution 
of the Q-ball equation, Eq. (|10|)) that the L-balls in this case will correspond to thick- 
walled L-balls, with radius ^ and charge L ^ (sm)~^, which, for typical value 



of s and m will not be much larger that 10^. This can be understood by noting that 
the effective tunnelling potential in this case is given by 

- f7(</') = ( Y + ^ - exp(-.0) j 02 . (13) 

Thus as ^ 0, this tends to ^{u^ — m^)^^. A tunnelling solution can only exist if 
w'^ < m?. In addition, the largest value of cannot be much larger than s~^. From 
this we can see that the largest charge of a L-ball of volume V has an upper bound 
L ^ ms~'^V . From the equation of motion it is easy to see that any solution of the 
L-ball equation will tend to zero once r ^ m~^. Thus L ^ (sm)"^ in general. Therefore 
the L-balls have a maximum charge and an essentially fixed radius in this case, with 
the scale set by the mass scale of the soft SUSY breaking terms. This is in contrast 
with the case of Q-balls which can have a thin-walled solution, for which the Q-ball 
charge is proportional to its volume for sufficiently large charge 0. 



3.3 u^dFdF and u^u^d^e^^ directions 



In these cases we consider the solution of the Q-ball equation for the potentials of 
Eq. d^) and Eq. (^. We find that the B-balls can have a thin- wall solution for suffi- 
ciently large charge, with a wall thickness of the order of uj^^. This can be understood 
as follows. For the case of the u'^df'd'^ potential, the equation for the B-ball is given by 



+ 



-uolcj) + m^0i^ log 




(14) 



where Uo is defined by 



2 2 



1 + 1 + l0£ 



(15) 



and 0c is the value of for which the right-hand side of Eq. ( p!4|) vanishes. If the initial 
value of is equal to 0c, then will remain constant for all r. However, if is slightly 
smaller than 0c, then will decrease slowly with increasing r until the change 50 is of 
the order of 0c, after which will quickly decrease to zero as r changes by 6r ^ 
(this gives the thickness of the wall of the B-ball) . The value of r at which begins to 
rapidly decrease, Tc, may be estimated by perturbing the B-ball equation around 0c. 
We find that 

(. 



(16) 



Thus, for sufficiently small 50(0), the B-ball can be made as large as we wish. This cor- 
responds to the thin-wall Q-ball property that the volume of the Q-ball is proportional 
to its charge. The value of ujo is determined by the solution of the B-ball equation for 



a given B. From solving numerically for the B-ball we find that uJq is approximately 
A\K\^/'^mo for \K\ ^ 0.1. Therefore, so long as the charge of the B-ball is sufficiently 
large that its radius in the thin-wall limit, Eq. (jlB]), is large compared with cu"^, we 
can use the thin-wall expressions. These are given by 




E 

Q 



and 

V ^ 



(17) 



(202t/(0„))i/2 ' V ; 

where 0o is the value of at the minimum of f/(0)/0^. If the charge of the B-ball is 
such that the radius from equation Eq. (|16D is not larger than uj~^, then the B-ball will 
be thick-walled. In this case, from numerical solutions of the B-ball equation, we find 
that as B decreases its radius stays more or less constant at roughly i/uo, but that 
the value of (j) inside the B-ball decreases. This is as we would expect, since the width 
of a bubble wall is purely determined by the mass of the associated scalar and we are 
effectively solving for a bubble here. The energy per unit charge in the thick-walled 
case is found numerically to be larger than in the thin-walled case. (In both cases it 
is not very much smaller than the free mass, in contrast with the case of gauge- 
mediated SUSY breaking where it is proportional to Q~^^^ 0). We may estimate the 
charge Qc at which a Q-ball becomes thin-walled by using the thin-wall expression for 
the Q-ball volume as a function of its charge, equation Eq. (|l8l), with the radius equal 
to {\K\^/'^mo)~^. For the d=6 u'^d^d^ potential Eq. (^, taking rrio ~ 100 GeV, we find 
that Qc ^ lO^^lir]-^/^, whilst for the d=4 u^'u^d'^e'' potential, Eq. (D, we find that 

The fact that the B-balls in these cases can be thin- walled is important, as it implies 
that there is no upper bound on the B-ball charge. This allows for the possibility of 
large charge and so relatively stable B-balls, which may have significant consequences 
for cosmology. 



4 Q-Ball Cosmology 

4.1 AD baryogenesis and Q-ball formation 

In general, we expect that the soft SUSY breaking mass squared term will receive 
an order correction, coming from non-minimal kinetic terms in the supergravity 
Kahler potential P]. Such terms would be expected in realistic unified theories such 
as string-type theories. The resulting potential for the complex AD field is then of the 



form 



f/(« « K - CsH-'M + + + . (19) 



where Ax = Ax^ + o,xH is the corrected A-term and cs and oa are typically of the 
order of 1. In this we are choosing the correction to have a negative sign in order 
to have condensate formation and AD baryogenesis. 

The minimal cosmological scenario for AD baryogenesis begins with a period of 
inflation, during which the Hubble parameter H has a value Hj ^ 10^^ GeV, in accor- 
dance with COBE perturbations [0. The infiaton subsequently oscillates coherently 
about the minimum of its potential until it completely decays and reheats the Uni- 
verse to a temperature Tr, which should be less than around 10^ GeV in order not 



to thermally regenerate too many gravitinos [0. The mass squared term of the AD 
scalar becomes positive at H ^ ms, at which time it begins to oscillate and the B 
or L asymmetry is formed by the influence of the phase dependent A-terms at this 
time. Throughout this period the Universe is matter-dominated by the energy of the 
coherently oscillating infiaton. For completeness, let us note that during infiaton os- 
cillation domination there will be a thermal background of particles, with temperature 
T Ri QA{MpHTlY/^, coming from infiaton decays 0. 

The basic picture for the formation of Q-balls within this scenario is as follows. As 
first observed by Kusenlco and Shaposhnikov |p, the AD condensate can be unstable 
with respect to growth of space-dependent perturbations of the condensate field. This 
is due to the fact that when the potential can lead to the formation of Q-balls, it 
must be increasing less rapidly that 0^ for some range of 0. If the condensate starts 
oscillating with an initial amplitude lying within this range of 0, then, on averaging 
over the coherent oscillations, it will behave as matter with a negative pressure. This 
negative pressure will cause any small perturbations to grow exponentially and go 
non-linear. It is assumed that the further collapse of the condensate in the non-linear 
regime will lead to the formation of Q-balls, once the non-zero charge (baryon or lepton 
asymmetry) of the condensate comes to stabilize the collapsing condensate field 0. 
This leaves the question of the origin of the initial "seed" perturbations. Since the 
AD field should be out of thermal equilibrium in order to avoid thermalizing the 
condensate before the baryon or lepton asymmetry can be created, the only source of 
initial perturbations in this scenario are quantum fluctuations of the AD field created 
during inflation. 

Thus, in order to discuss the length scale and charge of the regions of condensate 
which first go non-linear and collapse to produce Q-balls, we must first find the spec- 
trum of perturbations of the condensate field at H ^ ms due to quantum fluctuations 
and then consider the growth of perturbations of the condensate once begins to 
oscillate coherently at H ^ ms. 



4.2 Seed perturbations from quantum fluctuations during in- 
flation 



Let A be the length scale of a perturbation at H ms, when the AD condensate 
begins to oscillate coherently. Let Ag be the corresponding length scale at the end of 
inflation. Since the Universe is matter dominated after inflation, we have 

Ae \msJ 

The corresponding length scale at AA^e e-folds before the end of inflation will be 

AAATe = Aee-^^« . (21) 

This length scale will leave the horizon during inflation when Xan^ ~ H^^. Quan- 
tum perturbations of the A-D field become classical on leaving the horizon, with the 
magnitude of the perturbation, treating as a massless scalar, being given by 

~ ^ ■ (22) 

Here 6(f) refers to a perturbation about the minimum of the potential with H > ms, 
when has a negative mass squared term of the order of H^. Thus the physical 
scalar will have a mass of the order of H. Although this is not really an effectively 
massless scalar, which would require a mass small compared with H, it is not very 
massive relative to H and so we expect the massless result to be roughly correct. 
Once the perturbation is larger than the horizon it will appear to be a homogeneous 
perturbation on sub-horizon scales and so will coherently oscillate about the minimum 
of its potential and red-shift hke matter, 50 oc a~^/^. This assumes that the mass is 
large compared with H, so that we can ignore the 3i70 damping term in the equation 
of motion. Again this is not strongly satisfied, but we expect that it will roughly give 
the correct behaviour. Thus at the end of inflation the spectrum of perturbations will 
be given by 50e ^ e-^^^^^^y^ Hi /{2Tr), where AA^e(A) refers to the number of e-folds 
before the end of inflation at which the perturbation of length scale A first leaves 
the horizon. Suppose the perturbation re-enters the horizon at a^- Once back inside 
the horizon, the perturbation will be space-dependent. Its subsequent evolution will 
depend on whether its energy density is dominated by the gradient or potential energy. 
The 5(f) equation of motion is 

50 + 3HS^ - V^6(f) = -V'{6(f)) , (23) 

where V{6(f)) ~ H'^6(f)'^. The gradient energy will dominate if ^ H^, where 6(f> oc 
gjk x fjj^jg -^j}} obviously be satisfied once the perturbation enters the horizon. Since, 



as the Universe expands, is proportional to while is proportional to a~^, 
in general the perturbations will be gradient energy dominated once they re-enter the 
horizon. Therefore 5(f) oc a^^ . Putting all this together, we find that the spectrum of 
the sub-horizon sized perturbations at if ^ uis due to quantum fiuctuations during 
inflation is given by 



1 



^ Z ^1727777 • (24) 



4.3 Negative pressure and Q-Ball formation 

The growth of these perturbations once the AD field begins to oscillate coherently can 
be understood by considering the case of B-balls and the squark potential for values of 
somewhat smaller than the value of at the minimum of the potential (this describes 
the potential of the oscillating field once H is less than ms), 

(l + A-log(-J))^S (25) 

where we typically expect \K\ ^ 0.01 — 0.1. For lA'l small compared with 1, this has 
the form U{(f)) oc 0^^^'^. As discussed in reference |T^, the equation of state of matter 
corresponding to a coherently oscillating scalar in a potential of the form 0" is given 
by 

2n 

p=(7-l)p ; 7 = ^7^. (26) 

n + 2 



Thus for a scalar oscillating in the potential given by Eq. ( pS]) the equation of state is 
given by [T^ 



P = f P , (27) 

which gives a negative pressure if i^T < 0. Perturbations in the field will grow according 
to 

4 = , (28) 

where the perturbation of the density of the AD field on the length scale A = 27r/|k| 
is proportional to = Sp^/p. With 6p oc 2050, this implies that 

S<Pk = —S<P . (29) 

Thus the quantum fiuctuations of the field on the scale |k| = 27r/A will grow expo- 
nentially with time according to 

k = o(pik exp I I '^J^ ) t\ , (30) 



where t — corresponds to the beginning of the coherent oscillations and (pi and 50^ 
are the values of the field and its perturbation respectively at ~ ms- In this we 
have neglected the expansion of the Universe for simplicity. This is justified as we need 
only establish the value of H at which the perturbations go non-linear in an expansion 
time H~^. Since 50j and 0j scale the same way due the expansion of the Universe, we 
can calculate their ratio at if ~ 777.5. 

How can B-balls form as a result of this negative pressure? The naive picture 
is that a perturbation on some scale A will go non-linear once t ^ rrig^ (which is 
the smallest time scale on which B-balls can form and on which we can average over 
coherent oscillations to obtain a negative pressure), causing the AD condensate to 
collapse into fragments of size A, trapping some charge B, which will then form into 
B-balls with a charge of the order of B (depending on the efficiency of B-ball formation 
in the non-linearly collapsing fragments). However, by considering only the effect of 
negative pressure we are neglecting the effect of the charge B. The charge of the B-ball 
serves to prevent the soliton from collapsing to a radius smaller than the B-ball radius. 
Therefore we do not expect perturbations to be able to grow on length scales smaller 
than the B-ball radius once charge is taken into account. We can use the negative 
pressure arguement once the length scale going non-linear at a time t is larger than 
the final B-ball radius. We would then expect the B-balls to form quite efficiently 
when these B-ball size perturbations go non-linear. The charge of the B-ball will be 
given by the charge contained in a volume of radius of the order of the B-ball radius 
at this time. 

The time at which a perturbation of scale A goes non-linear is given by 




(31) 
where 

• <^^' 

We find that ak ~ 34 (44) for the d=4 (d=6) directions. In practice, the B-balls will 
typically turn out to be thick- walled, with radius Vo ~ {\K\^^'^ms)~^- Perturbations 
on this scale, which have the largest possible growth in a time H~^, will go non-linear 
at 

^^^T-' (33) 
\K\ms 

corresponding to if 0.1|ir|7775. In order to find the charge of the resulting B- 
balls we must find the charge density of the Universe at this time. We will assume 
that the charge asymmetry corresponds to the presently observed baryon asymmetry, 
rjB ~ 10~^°. Then the charge asymmetry at a given value of H during the infiaton 



oscillation dominated era prior to reheating is given by 

Thus we find that the charge of the region with radius ^ {\K\^^'^ms)~^ is given by 

Thus, with \K\ ^ 0.01, we see that B-balls of charge larger than 10^'* are likely to be 
formed, depending on the reheating temperature and the efficiency with which this 
charge is trapped within the B-balls. 

For the case of L-balls, we have seen that L-balls have a maximum charge of the 
order of 10^ and field strength of the order of 1 TeV, which is much smaller than the 
initial amplitude of the d=4 AD field at if ~ ms- This suggests that L-balls cannot be 
formed by the collapse of an unstable condensate at if ~ and that AD baryogenesis 
along the H^L direction will be essentially unaltered from the conventional scenario. 



4.4 Cosmological consequences of Q-balls in the MSSM 

Given that B-balls with a large charge will naturally form in the AD scenario along the 
squark directions, what might be their consequences for cosmology? This is crucially 
dependent upon the temperature at which the B-balls decay. If they decay after the 
electroweak phase transition has occured at Teui, then we will have the possibility of a 
new version of AD baryogenesis, in which the baryon asymmetry originates from the 
decay of the B-balls rather than the decay of the condensate. For example, if we were 
to consider L violating interactions which are in thermal equilibrium together with 
anomalous B+L violation at T ^ Tg^ (as may be expected in extensions of the MSSM 
with light Majorana neutrinos via the see-saw mechanism), then any B asymmetry 
coming from thermalization or decay of the AD condensate would be subsequently 
washed out [Q. However, the B contained in the B balls will not be erased, at 
least for some range of reheating temperatures, due to the fact that the large field 
inside the B-ball prevents thermal particles from penetrating beyond the surface of 
the B-ball, thus suppressing the thermalization rate To see this, we will give a 
conservative upper bound on the reheating temperature from requiring that the B- 
balls are not thermalized. Thermal particles will penetrate into the B-ball down to 
a stopping radius r^t, corresponding to g(j)st ~ T, where (p^t = (pi^st) and where the 
thermal particles gain a mass g(l) from interacting with the B-ball. Let us also assume 
that the thermal particles refiect from the "hard" region of the B-ball corresponding 
to r < Tst- Then there will be a fiux of thermal particles fiowing in and out of the 
"soft" part of the B-ball at r > r^t. The most effective way possible for this flux to 



remove charge from the B-ball would be via B absorbing inverse decay processes with 
a maximum possible rate Tinvd = kdT, where ~ 10^^ for strong interactions [^. 
The true inverse decay rate can be significantly smaller, depending on the details of 
the model. For large enough r, the thick-walled B-ball will be roughly described by 
a Gaussian profile, 0(r) ^ 0oe~^^/^°. Therefore the stopping radius will be given by 
Tst ~ rolog^^'^{(j)o/ 4>st)- Thus, since the charge density at radius r is ^^^^(r), which 
will be approximately constant over a change in radius 5r ^ rl/{2rst), we find that 
the maximum rate of loss of charge by the B-ball via inverse decays is approximately 
given by 

where ^ {\K\^^'^ms)~^ is the radius within which most of the Q-ball charge is 
concentrated. The condition for the Q-ball to survive thermalization is then that 
B/B ^ H. This gives an upper bound on the reheating temperature 

T„ ; ^^^^^ . (37) 

For thick-walled B-balls, the value of the scalar field at the centre of the B-ball, (po, is 
given by (po ~ 0.5\K\^^^mB^^'^, where B is given by Eq. (^). Thus, with g ^ 1 and 
kd ~ 10^^ for strong interactions and with m ~ lOOGel^, we obtain an upper bound 
Tr^ 5 X 10^\K\GeV . Thus, although a non-trivial upper bound may exist, even with 
relatively conservative assumptions there is still a range of reheating temperatures for 
which the B-balls can survive thermalization. Clearly a much more detailed analysis 
of thermalization is required to give the true upper bound on Tr, which could easily 
be weaker than our conservative estimate. (We have checked that destruction of the 
B-ball by direct collisions of thermal particles with the hard region is less important 
than thermalization of the soft region PB| ). 

From this we may conclude that, at least for some range of Tr, a B asymmetry can 
be preserved in the presence of rapid L violating interactions so long as the B-balls 
decay at < Tew Another interesting possibility is connected with the d=4 B-L 
conserving directions uVd'^e'^ and QQQL. These cannot produce a B asymmetry in 
conventional AD baryogenesis, since the condensate will be thermalized at tempera- 
tures large compared with Tg^ and so anomalous B+L violation will wash out the B 
asymmetry. However, if the Q-balls, which in this case are charged under U{1)b+l, 
decay after the electroweak phase transition, then a non-zero B asymmetry will result, 
together with an equal L asymmetry. In addition, even in those cases where conven- 
tional AD baryogenesis can, in principle, account for the observed baryon asymmetry, 
the B-ball decay mechanism might replace the conventional mechanism if, after the 
collapse of the AD condensate, most of the asymmetry were trapped in the B-balls. 



The Q-ball decay rate to hght fermions is proportional to the area of the Q-ball. It 
has been estimated to have an upper bound, which is hkely to be saturated for Q-balls 
with 0o much larger than nio, given by H] 



dQ uj'^A 



dt 1927r2 



(38) 



where A is the area of the Q-ball and where uj ~ ruo for small compared with 
1. The decay rate will depend on whether the Q-ball is thin or thick- walled. For 
thin-walled Q-balls the area of the Q-ball is related to the charge by 

^^M^. (3.) 
where k = (20^[/(0o))^/^. Thus the lifetime of the Q-ball in this case is given by 

For the thick-walled case, the area of the Q-ball is independent of its charge, being 
fixed by its radius Vo ~ {\K\^^'^mo)~^. The Q-ball lifetime in this case is then given by 

r - (41) 

For charges less than of the order of 10^^|i^|~^/'* for the u^d^d^ direction and 10^^|i^'|~^ 
for the u'^u'^d'^e'^ direction, the Q-balls will be thick-walled. The temperature at which 
the Q-balls decay is given by, assuming radiation domination, 

/ 1 \ .1/2 

T.-^) = . (42) 



1/2 



where kx = ( M . Thus, assuming ^ m-o ~ 100 GeV, thick- walled Q-balls 
will decay at 

' - ^"=^l°"^"%eV. (43) 



|ir|i/2 VlOO GeV/ \Q ) 

using /ct ^ 17. Thus we see that the Q-balls will decay at a temperature less than 100 
GeV if Q ^ 2 X lO^^I-f^l"^. For thin-walled Q-balls, the right-hand side of equation 
Eq. (|43|) has an additional factor {Q/Qc)^^'^, where Qc is the value of Q at which 
the thin-wall limit becomes valid. Q-ball decay is summarized in Fig. 2. Since we 
have shown that it is likely that the B-balls will have a charge greater than 10^^ in the 
MSSM, we see that it is quite natural for MSSM B-balls to decay after the electroweak 
phase transition has occured. Thus it is possible to generate the observed B asymmetry 
via the decays of B-balls occuring after the electroweak phase transition. This is true 
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Figure 2: Q-ball decay temperature T vs. the charge Q. The regions where L-balls 
and B-balls exist are also indicated. 

even in those cases with rapid L violating interactions or B-L conservation, where the 
conventional AD mechanism fails. 

The decay temperature depends on B, which in turn depends on the reheating 
temperature Tr. Nucleosynthesis requires that 1 MeV, which imposes an upper 
limit on the possible charge of the B balls, B ^ 10^^. This in turn imposes a lower 
bound on T^. For the d=4 case we find that ^ 50 MeV, whilst for the d=6 case we 
find Trf > 50 keV (see Fig. 2) . In principle, B-balls could decay at temperatures as low 
as that of the quark-hadron phase transition or less. This could have interesting con- 
sequences for the spatial distribution of baryons and inhomogeneous nucleosynthesis 

II- 

L-balls, which have Q ~ 1000, would decay at Td ~ 10^ GeV even if they could 
be formed primordially. It is therefore unlikely that such L-balls could have any 
cosmological consequences. However, it is possible that L-balls, which have a field 
strength of the order of 1 TeV or less, could play a role in the physics of the electroweak 
phase transition. 



5 Conclusions 

We have shown that it is quite natural for Q-balls to exist in the MSSM without requir- 
ing any special choice of SUSY breaking mechanism. For the case of B-balls associated 
with D-fiat directions involving squark fields, these are usually sufficiently stable to 
decay after the electroweak phase transition has occured. As a result, the Affleck-Dine 



mechanism along these direction will have to be revised, as the baryon asymmetry will 
come from condensate collapse and B-ball decay rather than simply from the decay of 
a coherently oscillating scalar condensate. As many extensions of the MSSM introduce 
additional L violating interactions, which, if they are sufficiently rapid to be in thermal 
equilibrium when anomalous B+L violation is in thermal equilibrium, rule out conven- 
tional AD baryogenesis but not B-ball decay baryogenesis, the potential importance 
of B-ball decay baryogenesis in the MSSM and its extensions is clear. We have also 
pointed out the existence of L-balls of small charge, characterized by a field strength 
not much larger than a TeV at most. Although such L-balls will be very short-lived, 
decaying at a temperature of around 10^ GeV even if produced primordially, and so 
are unlikely to significantly alter the conventional Affleck-Dine mechanism, they might 
nevertheless be produced during the electroweak phase transition. 

There are several points that this discussion of Q-balls in the MSSM has not ad- 
dressed. The details of the formation of B-balls from the non-linear collapse of the 
condensate and the efficiency which which charge is trapped in B-balls needs to be 
clarified, in order to obtain an accurate estimate of the B-ball charge, which deter- 
mines the temperature at which the B-balls decay. The physics of the thermalization 
of B-balls should also should be analysed in detail, in order to obtain constraints on 
the parameters of the models, in particular the reheating temperature. Given that 
the B-balls are decaying after the electroweak phase transition has occured, we should 
also consider what are the possible consequences of B-ball decay for other issues in 
cosmology. If the B-ball has a sufficiently large charge, corresponding to a low enough 
reheating temperature after inflation, then it might decay at a temperature close to 
that of the quark-hadron phase transition or nucleosynthesis. This could have con- 
sequences for inhomogeneous nucleosynthesis models. B-ball decay might also have 
consequences for dark matter in SUSY models, given that the B-balls are made of 
squarks which will produce neutralinos when they decay. For the case of L-balls, it 
would be interesting to consider the possible consequences of L-ball production during 
the electroweak phase transition. These issues will be the subject of future work. 
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